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The calculation of one-electron density of state values from the coefficient y of the term of the low 
temperature specific heat linear in temperature is complicated by many-body effects. In particular, the 
electron-phonon interaction may enhance the measured y as much as twofold. The enhancement factor 
can be evaluated in the case of superconducting metals and alloys. In the presence of magnetic mo- 
ments, additional complications arise. A magnetic contribution to the measured y was identified in the 
case of dilute alloys and also of concentrated alloys where parasitic antiferromagnetism is superim- 
posed on an over-all ferromagnetic order. No method has as yet been devised to evaluate this magnetic 
part of y. The separation of the temperature-linear term of the specific heat may itself be complicated 
by the appearance of a specific heat anomaly due to magnetic clusters in superparamagnetic or weakly 
ferromagnetic alloys. 
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1. Introduction 

In the Sommerfeld-Bethe theory of metals the elec- 
tronic specific heat at low temperatures is linear in tem- 
perature in first order approximation. The lattice 
specific heat in the low temperature approximation is 
proportional to T 3 so that, in the absence of other con- 
tributions, the total specific heat 



C = yT + j$T* 



(1) 



If C is known as a function of T, the two terms can be 
separated by making use of the linear variation of CjT 
with T 2 and by extrapolating to T=0. The intercept of 
the extrapolated line with the ordinate axis gives the 
temperature coefficient of the electronic specific heat 
y. In the simplest case, y is proportional to the elec- 
tronic density of states at the Fermi surface, N(Ef): 

y = (H3)7T 2 ak 2 N(E F ) (2) 

where k is the Boltzman constant and a is a numerical 
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factor determined by the units used for y, N(Ef) and k. 
Unfortunately, in a very large majority of cases, the 
simple procedure just described cannot be used, or at 
least it does not give reliable results. Many-body effects 
and, in some alloys, magnetic effects may make the 
determination of N(E F ) from low temperature specific 
heat data more complicated than implied by eqs (1) and 
(2), or even impossible at the present state of the art. 

2. Many-body Effects 

In recent years it has become known that many-body 
effects, in particular the electron-phonon interaction, 
require renormalization of the effective mass of the 
electrons at the Fermi surface. This increases the mea- 
sured electronic specific heat coefficient over the one- 
electron "band structure" value by the enhancement 
factor (1 + \). For Na, Al and Pb, it was possible to 
determine the value of this factor [ 1] , by comparing the 
"band structure electronic specific heat," calculated 
from the known band structure and the topography of 
the Fermi surface, with the measured electronic 
specific heat. These values: 1.25, 1.45 and 2.00, respec- 
tively, were found to agree quite well with the enhance- 
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ment factors calculated from band structure, Fermi 
surface topography and phonon dispersion curves, on 
the basis of the electron-phonon interaction [1] . Unfor- 
tunately, for most other metals calculations of this sort 
cannot be made at present since at least some of the 
required data are not yet available. For superconduct- 
ing metals the electron-phonon coupling constant k has 
been recently calculated by McMillan [2] , using the 
following equation which he derived from the strong 
coupling theory: 



T c = 







1.45 



exp- 



1.04(1 + A.) 
\-/x*(l + 0.62X) 



(3) 



where Tc is the superconducting transition temperature 
and 6 is the Debye temperature. The electron-electron 
interaction constant /jl* was assumed to have a value of 
0.13 for all transition metals. The values of X calculated 
by McMillan [2] for superconducting metals are given 
in table I. 

Table I. The Electron-Phonon Interaction Coefficient k and 
"Band Structure" Density of States N(E F ) for Superconducting 
Metals. 2 



Element 



Be. 
AL. 

Zn. 
Ga. 
Cd. 
In.. 
Sn. 
Hg. 
TL. 
Pb. 
TL. 
V.. 
Zr.. 
Nb. 
Mo 
Ru. 
Hf. 
Ta. 
W.. 
Re. 
Os. 
Ir.. 



T c 
°K 


e 

°K 


A 


N(E F ) 
eV _1 atom 


.026 


1390 


.23 


.032 


1.16 


428 


.38 


.208 


.85 


309 


.38 


.098 


1.08 


325 


.40 


.091 


.52 


209 


.38 


.106 


3.40 


112 


.69 


.212 


3.72 


200 


.60 


.238 


4.16 


72 


1.00 


.146 


2.38 


79 


.71 


.182 


7.19 


105 


1.12 


.276 


.39 


425 


.38 


.51 


5.30 


399 


.60 


1.31 


.55 


290 


.41 


.42 


9.22 


277 


.82 


.91 


.92 


460 


.41 


.28 


.49 


550 


.38 


.46 


.09 


252 


.34 


.34 


4.48 


258 


.65 


.77 


.012 


390, 


.28 


.15 


1.69 


415 


.46 


.33 


.65 


500 


.39 


.35 


.14 


420 


.34 


.51 



It is now clear that the enhancement factor (1 4- k) of 
the electronic specific heat due to the electron-phonon 
interaction can be as high as 2, or more. This interac- 
tion affects only the electrons whose kinetic energy is 
close to the Fermi energy. The density of states at lower 
levels, that is for most of the electrons in the metallic 
band, may be assumed to correspond to the one-elec- 



tron "band structure" situation. Hence, the lower, 
"band structure density of states" values must be used 
in determining the band width, for instance, rather than 
the density of states enhanced by electron-phonon in- 
teraction, as obtained from low temperature specific 
heat measurements. Using the electron-phonon 
coupling constant X, for instance the values given in 
table I, the "band structure density of states" at the 
Fermi level N(Ef) can be calculated from the experi- 
mentally determined value of the low temperature 
specific heat coefficient y' as follows: 

3 



N(E F y- 



'rt 2 «(l + A) 



y 



(4) 



For most of the nonsuperconducting metals and alloys 
the value of k is at present unknown and, as a result, 
the "band structure density of states" cannot be calcu- 
lated from the low temperature specific heat. 

As seen in figure 1, the experimental electronic 
specific heat coefficient y' for the b.c.c. 3d-transition 
metals and their alloys as a function of electron concen- 
tration [3] shows prominent maxima and minima in the 
range of e/a from 4 to 9. Since in the region of the 
minima and of the second maximum the alloys are not 
superconducting, the "band structure density of 
states" cannot be calculated at present. Thus, the in- 
teresting question whether the prominent features of 
these curves are due to changes in the electron-phonon 
enhancement factor upon alloying, or indeed these fea- 
tures are characteristic of the electronic band structure 
of the transition metals concerned, cannot be answered 
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Figure 1. Coefficient y' or y" of the low temperature specific heat 

term linear in temperature vs electron concentration e/a for b.c.c. 

alloys of 3d transition metals [3] . Points marked by filled squares 

represent data for close-packed structures. 
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with certainty. However, the work of McMillan [2] al- Tg i^ 
lows the conclusion that the electron-phonon coupling 
constant (and, thus, the enhancement factor) depends 
primarily on the phonon frequencies, rather than on the 
electronic properties. Since the elastic constants and, 
therefore, the phonon frequencies are not known to 
undergo drastic changes with the composition in such 
solid solution alloys composed of metals near one 
another in the same row of the periodic table, it may be 
concluded with a reasonable degree of probability that 
the prominent features mentioned of the y' versus e/a 
curve of figure 1 are in fact resulting from correspond- 
ing variations in the "band structure density of states," 
even though the relative magnitude of the various 
minima and maxima may be appreciably altered by the 
gradual changes in the coupling constant with composi- 
tion. 
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3. Magnetic Effects 



Considerable difficulties are often encountered in 



T 2 in DEG 2 

FIGURE 2. (C-A)/T vs T 2 (where C-A is low temperature specific heat 

less magnetic cluster contribution, see eq (5)) for alloy NIo.ah Cuo.sa 

cooled without a magnetic field (top graph), cooled in 14 kOe field 

from 300 to 4.2 K with the field turned off during the measurements 

(graph m) and fie Id-cooled with the field on during measurements 

(graph m — m) [8]. 

The occurrence in the same alloys of "exchange 
anisotropy" and of a magnetic contribution to the low 
temperature specific heat term linear in temperature 
suggests that these two phenomena may be associated 
with the same structural condition. This expectation is 
further supported by the fact that y m is normally also af- 
fected by field cooling. According to Kouvel's highly 



determining the value of y' for solid solution alloys of successful model [7] , the structural condition responsi- 
ferromagnetic with antiferromagnetic or nonmagnetic 
metals. For instance, it was found [4] for the random 
solid solution alloys Mn-Ni that, in addition to the elec- 
tronic specific heat coefficient y', the measured coeffi- 
cient y" of the term of the low temperature specific he ^ 
linear in temperature includes also a magnetic con- 
tribution y m : 



y /,=< y' + 7m. 



(5) 



The alloy MnNi ;} can be ordered by thermal treat- 
ment and, in the well-ordered condition, the coefficient 
of the linear term of the low temperature specific heat 



ble for "exchange anisotropy" is a spatially in- 
homogeneous magnetic state, e.g., the superposition of 
local "parasitic antiferromagnetism" on net overall fer- 
romagnetism. Overhauser [9] and Marshall [10] con- 
nected the magnetic contribution to the linear term of 
the low temperature specific heat with the location of 
a sufficient number of spins in a near-zero field. In 
Overhauser's theory this condition arises at the nodes 
of the static spin density waves of an antiferromagnet. 
Marshall pointed out that the required condition may 
arise in dilute spin systems, where the average distance 
between neighboring spins is sufficiently large, so as to 
make the interactions weak, as in dilute Cu-Mn alloys. 



is approximately half that for the disordered alloy of the The alloys considered above are neither antiferromag- 
same composition. This lower value is substantially netic nor dilute. However, because of the peculiar, 
free of the magnetic contribution y m and it may be con- complicated spin arrangement, resulting from the su- 
sidered as approximately equal to the real experimental perposition of local parasitic antiferromagnetism on net 
electronic specific heat y' of the alloy. On the other overall ferromagnetism, it may be expected that many 
hand, the larger y" value for the disordered alloy in- spins are located in small regions where ferromagnetic 
eludes y m . Similar magnetic contributions to y" were and antiferromagnetic exchange interactions nearly 
identified in a number ofotherf.ee. solid solution alloy cancel each other locally, so that the average field in 
systems [4] and in b.c.c. Fe-Al alloys [5]. It is signifi- such regions is near zero [4]. In accordance with this 
cant that in the same alloy systems, and at similar com- "local field-cancellation" model, the effect of field cool- 
positions, magnetic measurements by Kouvel [6,7] de- ing on y m may come about if the application of an exter- 
tected the appearance of an asymmetrical hysteresis nal magnetic field during cooling through the Curie 
loop after cooling in a magnetic field through the Curie temperature increases or decreases the number of 
temperature ("exchange anisotropy"). In addition to this spins located in near-zero field. Both increase and 
effect of field cooling on the magnetic properties, in decrease [(MnNi ;i [4], Nio.4sCuo.52 [8])] were in fact 
several instances an effect of field cooling on y m was observed. It is easy to visualize that the change in y m as 
also detected [4,8] , figure 2. a result of field cooling may also happen to be negligibly 
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small, even though the value of the magnetic contribu- 
tion y m itself is large. Thus, while the occurrence of a 
measurable effect of field cooling on the temperature- 
linear term of the low temperature specific heat may be 
considered as a proof for the existence of a magnetic 
contribution to this term, the absence of such an effect 
does not prove that y m is zero, or that it arises through 
a mechanism different from the "local field cancella- 
tion." 

If the experimentally determined coefficient of the 
temperature-linear term of the low temperature 
specific heat includes a magnetic contribution, it is at 
present not possible to derive from such a y" value the 
"experimental electronic specific heat coefficient" y\ 
which is free from y m . This is well illustrated by the Ni- 
Cu alloys, for which the coefficient of the linear term 
has a maximum around the composition Nio.4sCuo.52 
[8]. A detailed study of the properties of these f.c.c. 
solid solutions at compositions in the vicinity of the 
maximum [11] shows that the experimental values of 
the coefficient do in fact include a magnetic contribu- 
tion. It is, therefore, not possible to tell whether the 
maximum is entirely due to y m , or whether y f itself has 
a maximum, which is merely increased by the addition 
of y m . A maximum in y' has been expected on theoreti- 
cal grounds [12] because of enhancement due to the 
electron-paramagnon interaction [13,14]. The theory 
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FIGURE 3. Coefficient y" of the temperature-linear low temperature 

specific heat term, Curie temperature TV, temperature independent 

specific heat term A and coefficient /3 of the T 3 term, obtained by least 

squares fitting to eq (5), for Cu\. x Ni x (solid lines) and for (Cu\. x Ni x )o.9 

Alo A (dashed lines) alloys vs x [11 ,22,24] . 




g 



30 

20 
A 
10 



-0.1 







*\ 


1 


ft 


/ 


* \ 

iv 


o> 
a> 


^R 


r 


\ 

\ 
• 


b ' 

— 0) 


A \ 


i 


\ 


S. 


i» 






f*% 


*^-o 


f 




- ) 


<t 


k* y 


\/ 








^•' 





0.2 



0.3 



0.4 



0.5 



0.6 



FIGURE 4. Low temperature specific heat and Curie temperature data, 

as in figure 3, for V\. x Fe x (solid lines) and (V\. x Fe x )o.9 AIqa (dashed 

lines), alloys [3 ,23,24]. 

would require the maximum of y' to occur at the critical 
composition, where ferro magnetism just begins to set 
in at K. Detailed study of the magnetic properties 
showed [11] that the critical composition is at approxi- 
mately 57 percent Cu, where y" already decreased to a 
value far below its maximum. Figure 3 gives y" for the 
Ni-Cu solid solutions, together with Curie temperature 
data, which define the critical composition. One may 
conclude that the maximum in y" is largely, or entirely, 
due to y m , rather than to y' . That the maximum in y m 
should occur in the weakly ferromagnetic region is en- 
tirely consistent with the "local field cancellation" 
model discussed above. It is quite likely that the max- 
imum in the coefficient of the linear term for the f.c.c. 
Rh-Ni solid solutions, which apparently also occurs off 
the critical composition on the ferromagnetic side 
[15,16] , is also due to a magnetic contribution y m . In 
fact, several solid solution alloy systems are now known 
to exhibit similar conditions. Figure 3 shows this for a 
series of Ni-Cu-Al ternary alloys with a constant Al-con- 
tent of 10 percent. It is seen that the Al addition shifts 
the critical composition to a higher Ni/Cu ratio, that the 
maximum in y" (which is here even higher than for the 
binary alloys) is also shifted, and that it again appears 
away from the critical composition, on the ferromag- 
netic side. Further examples are given in figure 4, 
which shows similar data for b.c.c. V-Fe binary and V- 
Fe-Al ternary solid solutions with a constant Al-content 
of 10 percent. 
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As yet there appears to be no experimental evidence 
for the theoretically predicted [14] peak in y' in alloys 
at the critical composition, resulting from the electron- 
paramagnon interaction. 

Figures 3 and 4 illustrate also a complication, which 
arises quite frequently in extracting the coefficient of 
the temperature-linear term from low temperature 
specific heat data for weakly ferromagnetic and almost 
ferromagnetic alloys. In many alloy systems in a certain 
region around the critical composition an anomaly in 
the measured low temperature specific heat is ob- 
served, so that C is no longer given by eq (1). This 
anomaly is conspicuously evident in the usual CjT 
versus T 2 graph. Instead of being a straight line, this 
graph becomes a curve, extending upward at low tem- 
peratures. It was shown by Schroeder [16] that, in 
such cases, eq (1) can be usually replaced by 



10" 



C=A+yT + ($T\ 



(6) 



Schroeder found that the addition of a temperature-in- 
dependent term A to the low temperature specific heat 
results from the thermal excitation of magnetic clusters 
present in many nearly ferromagnetic and weakly fer- 
romagnetic alloys. The presence of magnetic clusters 
in Ni-Cu alloys in the composition range 50-56 percent 
Cu has been recently beautifully documented by Hicks, 
Rainford, Kouvel, Low and Comley [17] by means of 
neutron magnetic diffuse scattering. In Schroeder's 
theory the magnetic clusters, which interact with a 
weak crystal field, are thermally excited and they con- 
tribute an Einstein specific heat, which is temperature- 
independent above the Einstein temperature. The tem- 
perature range of 1.4 to 4.2 K, frequently used in low 
temperature specific heat measurements, appears to be 
above the Einstein temperature in most such systems. 
Investigations by Scurlock [18] show a decrease at 
lower temperatures of the anomalous specific heat from 
its constant value above 1.4 K, suggesting that the Ein- 
stein temperature is near 1.4 K. Since A includes an 
equal contribution of k for each cluster, regardless of 
the cluster moment [ 16] , the low temperature specific 
heat data, from which the value of A can be extracted 
by least squares fitting to eq (6), give reliable informa- 
tion as to the number of thermally excited magnetic 
clusters. The correlation of the number of clusters for 
Ni-Cu alloys by specific heat measurements [11] with 
the number of clusters derived from neutron scattering 
[ 17] is given in figure 5. It is seen that the low tempera- 
ture specific heat data for the superparamagnetic alloys 
(Cu-content larger than 57%) are quite consistent with 
the neutron scattering data. For the weakly ferromag- 
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Figure 5. Concentration of magnetic clusters in Ni-Cu alloys vs 
composition [22~\ . Filled circles represent cluster concentrations 
calculated from specific heat data. Empty circles give magnetic 

cluster concent rations derived from neutron magnetic diffuse 
scattering [17]. 



netic alloys the number of thermally excitable clusters 
rapidly decreases with increasing Curie temperature 
(decreasing Cu-content). In this composition range 
most of the magnetic clusters interact with one another 
and become part of the ferromagnetic system. As seen 
in figures 3 and 4 the maximum of A corresponds well 
with the critical composition for all four alloy systems 
considered. Figures 3 and 4 also show the anomalous 
behavior of the parameter /3, obtained by least squares 
fitting to eq (6). The anomalous variation of ft with com- 
position occurs in all four alloy systems in the vicinity 
of the critical composition, and it is clearly magnetic in 
origin [5]. 

4. Nuclear Specific Heat Effects 

It was shown by Marshall [19] that the hyperfine in- 
teraction between the dipole moments associated with 
certain nuclides and the effective field H eJ y at these 
nuclei, resulting from electronic spin moments, gives 
rise to a contribution to the low temperature specific 
heat. This contribution decreases rapidly with increas- 
ing temperature; in the 1.4 to 4.2 K range it is propor- 
tional to T~ 2 ("high temperature" approximation). Con- 
sequently, the nuclear magnetic specific heat term can 
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378-324 O - 70 - 10 



be separated easily from the electronic specific heat 
term yT in that temperature range. The nuclear 
quadrupole specific heat recently reported by Phillips 
[20] and by Martin [21] is also proportional to T' 2 and, 
thus, poses no problem in determining the electronic 
specific heat. 
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